Particle density distributions in Fermi gas superfluids: Differences between one and two channel 

models 
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We show how to describe the T / behavior associated with the usual BCS- Bose Einstein condensation 
(BEC) crossover ground state. We confine our attention here to the BEC and near-BEC regime where analytical 
calculations are possible. At finite T, non-condensed fermion pairs must be included, although they have been 
generally ignored in the literature. Within this BEC regime we compute the equations of state for the one 
and two channel models; these two cases correspond to whether Feshbach resonance effects are omitted or 
included. Differences between these two cases can be traced to differences between the nature of a Cooper pair 
and bosonic condensate. Our results are also compared with the Gross Pitaevskii equations of state for true 
bosons. Differences found here are associated with the underlying fermionic character of the system. Finally, 
the particle density distribution functions for a trap containing superfluid fermionic atoms are computed using a 
Thomas-Fermi approach. The one and two channel behavior is found to be very different; we find a narrowing 
of the density profile as a result of Feshbach resonance effects. Importantly, we infer that the ratio between 
bosonic and fermionic scattering lengths depends on the magnetic detuning and is generally smaller than 2. 
Future experiments will be required to determine to what extent this ratio varies with magnetic fields. 



I. INTRODUCTION 

The recent observations ClIlllBlllSBIlof 
Bose-Einstein condensation (BEC) of molecules formed from 
fermionic atoms are extremely exciting. Because of a 
Feshbach resonance, in which fermions couple to molecu- 
lar bosons, it is possible 0, Uoll . via application of mag- 
netic fields, to obtain an arbitrarily strong attraction between 
fermions, and to probe the crossover II ll ll 211 from BEC to 
BCS. The resultant superfluidity of pre-formed pairs has a nat- 
ural counterpart in some theories I13l fl4l \19L Hal of high T c 
superconductors. Indeed, a striking, and important feature of 
the cuprate superconductors is their pronounced precursor su- 
perconductivity, or "pseudogap" effects, the origin of which 
is still under active debate. 

A second, very important motivation for these experi- 
ments is based on the theoretical observation that a BCS-like 
ground state wavefunction is capable 11711 of describing both 
fermionic and bosonic -based superconductors, provided that 
the chemical potential of the fermions, /x, is determined self 
consistently. Given the vast success of weak coupling or BCS 
theory, it is extremely important to formulate this extended 
theory at all T, and confront it with controlled experiments. 

In this paper we explore the implications of this specific 
ground state wavefunction 0111 in detail to address both 
T = and T = T c , with particular emphasis on the "near- 
BEC" regime. Our goals are (i) to emphasize how to in- 
clude finite temperature effects in a manner consistent with 
the ground state; this necessitates the introduction of non- 
condensed fermion pairs which lead to "pseudogap" effects 
11911 . (ii) to provide analytic calculations and insights by 
working in a (near-BEC) regime where calculations are more 
tractable, (iii) to discuss in some detail the differences be- 
tween the "one channel" and "two channel" models for the 
Feshbach resonance, and (iv) to compare with well established 
theories of weakly interacting Bose superfluid 1 20], as well as 
(v) with the measured density distributions in a trap. We begin 



with the homogeneous case and then consider the trap con- 
figuration; from this we infer the ratio between the effective 
bosonic and fermionic scattering lengths, which is found to 
be strongly dependent on the magnetic detuning v , and gen- 
erally less than 2, as a consequence of many-body combined 
with Feshbach resonance effects. 



II. T MATRIX FORMALISM: BEYOND NOZIERES AND 
SCHMITT-RINK 

The Hamiltonian used in the cold atom and high T c 
crossover studies consists, in its most general form, of two 
types of interaction effects: those associated with the direct 
interaction between fermions parameterized by U, and those 
associated with "fermion-boson" interactions, whose strength 
is governed by g. 
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Here the fermion and boson kinetic energies are given by 
e k = h 2 k 2 /2m and e™ b = h 2 q 2 /2M, respectively, and v 
is an important parameter which represents the "detuning". 
The variational ground state which we will consider here is a 
product of both fermionic and bosonic contributions 



* = *0 ® 



(2) 



where the normalized fermionic wave function is the standard 
crossover state fT7LIT3l 
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and the normalized molecular or Feshbach boson contribution 
is represented by 



^ = e -A 2 /2+A6t| ) 



(4) 



The variational parameters are, thus, Uk, Vk and A. Applying 
standard variational techniques on the ground state wavefunc- 
tion leads to Eqs. @ and E3> and the T = limit of Eq. d26l . 
along with the result that A = <f> m = (6 q=n ). 

Whether both forms of interactions (U and g) in the Hamil- 
tonian are needed in either system is still under debate. The 
bosons (6 k ) of the cold atom problem |9, 10] will be referred 
to as Feshbach bosons. These represent a separate species, not 
to be confused with the fermion pair (a k at k ) operators. Thus 
we call this a "two channel" model. 

For the sake of clarity we begin by ignoring the Fesh- 
bach resonance-induced interactions. The variationally in- 
duced constraints on the fermionic degrees of freedom fPUl 
are given by 



= 1 



(5) 



k u k> 



where v k 



with the fermion density n = 
[1 — (ek — /i) /Ek) /2. We define the quasiparticle dispersion, 
Ek, below. 

One can anticipate that the excited states of this system are 
associated with their natural extension fl5l Q(| fl^l following 
BCS theory. For T < T c , the gap equation and number equa- 
tions now involve Fermi functions f(Ek) and are given by 
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w here /i is now evaluated at finite T and Ek = 
V ( e k ^ m) 2 + A 2 . Here we assume a contact potential where 
a s is the isolated two body scattering length, which is simply 
related to U, (see, for example, Eq. JA9» . 

Two important observations need to be made. Here A is 
associated with the excitation gap for fermions. This exci- 
tation gap reflects both condensed and non-condensed pairs 
and is not the same as the order parameter A sc , except at 
T = 0. Eqs. (|6jl and are widely cited in the literature, 
but with an important difference. Throughout the literature A 
is taken to be the order parameter. By contrast, to determine 
A sc here we will first need to characterize the contribution of 
non-condensed pairs to A. 

Our second important observation is that for T < T c , the 
Fermi factors in the above equations are essentially negligi- 
ble 1 2 1 ] . It follows that both A(T) and fi(T) are temperature 
independent in this near-BEC regime. Indeed, this is consis- 
tent with the physical picture of well established, pre-formed 
pairs in the BEC limit, so that the fermionic energy scales are 
unaffected by T below T c . 

The simple physics in Eq. l|6} may be schematically rep- 
resented by plots of A versus temperature. Figure^contrasts 




BEC 



FIG. 1: Contrasting temperature dependences of A in the BCS and 
BEC regimes. Similarly, in the BEC regime n is a constant, so that 
all fermionic energy scales are T independent, as expected. 



the behavior in the weak coupling or BCS and strong coupling 
or BEC regimes. In the BCS limit A(T) follows the behavior 
of the order parameter, whereas in the BEC regime, pairs are 
pre-formed and there is no temperature dependence in A(T) 
on the scale of T c . We now extend these qualitative observa- 
tions to a more quantitative level. 



A. Extending Conventional Crossover Theory to T ^ 0: BEC 
Limit Without Feshbach Bosons 



The self consistent equations in the BEC limit for general 
temperature T can then be written as 
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T < T c . (9) 



It follows from the above equations that, just as in the T = 
limit I22II23II . we have for general T <T C , 



Upairs 2 -^oA 

where the coefficient of proportionality 

Zo 



m 2 a s 
8tt?i 4 



(10) 
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We arrive at an important physical interpretation. Even 
though A or n pa i rs is a constant in T, this constant must be 
the sum of two temperature dependent terms. Just as in the 
usual theory of BEC these two contributions correspond to 
condensed and non-condensed components 



^ condensed 



f r r\ 1 ^non— condensed /rri\ /ion 
( T )+ n p al rs ( T )> (I 2 ) 



so that we may decompose the excitation gap into two contri- 
butions 



A 2 = A 2 C (T) + Ap (T) , 



(13) 



where A SC (T) corresponds to condensed and A pg (T) to the 
non-condensed (or pseudo) gap component. Each of these are 
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proportional to the respective number of condensed and non- 
condensed pairs with proportionality constant Z . Just as in 
BEC, at T c , 



^non— condensed /rri \ " \ A i/r\ rri \ 

" ... . {Tc) = 2 = Z^ b ( il i> Tc ) ■ 

q 



(14) 



where b{x) is the usual Bose-Einstein function and fl q is the 
dispersion of the non-condensed pairs, which will be self con- 
sistently determined below. Thus 

A 2 (T C ) = A 2 pg {T c ) = Z^J2W q ,T c ) = ~Z^ . (15) 



We may deduce directly from Eq. (I15> that A 2 g — 
— J^q t{Q)< ^ we P resume that below T c , the non-condensed 
pairs have propagator 



t{Q) = 



Zn 



(16) 



This lead to a key question: how can one deduce the contri- 
bution from non-condensed pairs? We now work backwards 
to infer the dispersion tt q for these pairs. A fundamental re- 
quirement on non-condensed pairs in equilibrium with a Bose 
condensate is that their effective chemical potential satisfy 



(17) 



Equation (|6) can be shown to be consistent with Eq. dl7> pro- 
vided that the propagator for non-condensed pairs is given by 



t{Q) 



u 



i + u x {Q) : 



where 



x{Q) = J2 Q { K ) G o{Q-K), 



(18) 



(19) 



K 



and G represents the fermionic Green's function which has a 
self energy E(-ff) = — A 2 Gq{— K). Here Go is the bare prop- 
agator. The details of this analysis are presented in Appendix 

m 

Another important point should be noted. This pair propa- 
gator or T-matrix differs from that first introduced by Nozieres 
and Schmitt-Rink ll24ll because here there is one dressed and 
one bare Green's function. In the approach of Ref. 24 both are 
taken as bare Green's functions. By contrast there are other 
schemes in the literature 1 25112a 12711 where both Green's func- 
tions are dressed. We end by noting that at small four-vector 
Q (and moderately strong coupling) we may expand Eq. Jl 8i 
to obtain 



t{Q) 



7 -i 



(20) 



Now we can deduce directly from Eq. Jl 8i that the dispersion 
of non-condensed pairs is of the form 



fl q = h 2 q 2 /2M* 



(21) 



In summary, this quadratic dispersion can be derived from the 
pair susceptibility x{Q)- in turn, the particular form for x(Q) 
shown in Eq. d!9l > is chosen in order to be consistent with 
Eqs. © and Q. In this sense the usual BEC constraint ex- 
pressed by Eq. (I17> is intimately connected to the BCS-like 
gap equation of Eq. (|SJl. The details of this analysis are shown 
in AppendixlAl 



B. Differences Between One and Two Channel Models 

We now extend this analysis to include Feshbach bosons 
For this situation we can write down an equation 
11911 equivalent to Eq. (|6j with the direct fermion interaction U 
replaced by U e f / = U + g 2 /(2/z — v) and effective scattering 
length a s — > a* — mil* / '{AttH 2 ). Here we define 



U* = U - 



{vq — 2/i) m 
where a* is dependent on fi. We thus have 

A(T) = -t/ e// ]TA(T)!- 2/(£;k) 
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where E k = ^{e^ - fi) 2 + A 2 (T), and again A is to be dis- 
tinguished from the order parameter 11911 . Alternatively one 
can rewrite Eq. (I23> as 

' 1 i-2f(E k y 

2e k 



Aiih 2 a* 



E 

k 



2£ k 



(25) 



In Eq. \2A\ n represents the number of fermions, but the 
number equation for the total number of particles involves 
both condensed and uncondensed bosons as well 



where n° 



n + 2n b + 2nl = n tot , (26) 

b 2 m is the number of molecular bosons in the 
condensate. The number of noncondensed molecular bosons 
is given by 

n b{T) = -Y< D (Q) ( 2? ) 
where the Bose propagator is 

1 



D{Q) = 



Sb(0) 



„ , + 2/i 
and we choose the self energy 1 19] 

X B (Q) = ~g 2 X {Q)/[l + U X {Q)} 



(28) 



(29) 



to be consistent with the Hugenholtz-Pines condition that 
bosons in equilibrium with a condensate must necessarily 
have zero chemical potential. This is equivalent to 



(J>boson{T) = 0, T <T C 



(30) 



where Hboson — 2/i — j/ + Eb(0). It follows after some simple 
algebra that Eq. d30i is then consistent with Eq. d23l >. 
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C. Parameter Choices and Renormalization Scheme 

We have introduced quantities (£/o, 30,^0) which charac- 
terize the two body scattering in vacuum and correspond to 
the scattering lengths a s or a*. How are these related |28] 
to the parameters (U, g, v) which enter into the Hamiltonian? 
To begin we ignore Feshbach effects. The way in which the 
two body interaction U enters to characterize the scattering 
(in vacuum) is different from the way in which it enters to 
characterize the N-body processes leading to superfluidity. In 
each case, however, one uses a T-matrix formulation to sum 
an appropriately selected but infinite series of terms in U. For 
the two-body problem in vacuum, there is a measurable char- 
acteristic of the scattering, the scattering length, a s , 



m 



1 



(31) 



which is related to U via the Lippmann-Schwinger equation 

1 X 1 



(32) 



In this way one may solve for the unknown U in terms of 
a s or Uq. The two-body T-matrix equation (Eg. (132^ can be 
rewritten as 



u = ru Q , r 



where we define the quantity U c as 



Uo 
U c 



- = -E- 



2e k 



(33) 



(34) 



Here U c is the critical value of the potential associated with the 
binding of a two particle state in vacuum. Specific evaluation 
of U c requires that there be a cut-off imposed on the above 
summation, associated with the range of the potential. 

Now, for the Feshbach problem, the calculation proceeds in 
a similar fashion and one finds |28] 

U = TU , g = Tg , v - v Q = -A , (35) 

V c 

Here we define the parameter v$ which is directly related to 
the difference in the applied magnetic field B and Bq 



vo = (B- B )A^°, 



(36) 



where AfjP is the difference in the magnetic moment of the 
two paired hyperfine states. To connect the various energy 
scales, which appear in the problem, typically IGauss w 
60.E F for 40 K. 



D. Characterizing the Condensate 

There is an important difference in the nature of the con- 
densate for the two cases, with and without Feshbach bosons. 



In the latter case there are two components to the condensate 
associated with Cooper pairs (A sc ) and condensed molecular 
bosons. Moreover, the Cooper condensate enters only into the 
gap equation, whereas the molecular Bose condensate enters 
also into the number equation. As long as the bosonic conden- 
sate n° is non-negligible, this difference leads to essentially 
different physics between the one and two channel problem, 
as will be demonstrated below. 

In this section we enumerate the different relationships be- 
tween the Cooper condensate and bosonic condensate terms. 
The order parameter associated with Eq. Q represents a lin- 
ear combination of both paired fermions (Coope r condensate) 
and condensed molecules. It is given by lTllfl2Tl 

A sc = A sc - g<p m (37) 

where, the boson order parameter <f> m = (6 q=0 ). We have 

A sc = -Z7^(a_ k; a kT ), (38) 



U 



therefore, A sc can be written as 

g 2 
(2/i - v) 

where we have used the fact l28ll that 

, gA sc 



^(a_ ki a kT ) 
k 



{v-2n)U 



(39) 



(40) 



The number of condensed Feshbach bosons which enters 
the number equation [Eq. (I26H is given by n° = <fi 2 n . Thus 
we have 



g 2 A 2 



g 2 A 2 



b [{v-2n)Uf [(2/i - u)U + g 2 } 2 
Using Eqs. j35l > we can also write 



„" = 1 _ Z» 

U* 



UoVA 2 



9l 



(41) 



(42) 



where U* is defined in Eq. \22\ . 



III. EQUATIONS OF STATE AT T = 

A. One Channel Model 

We now rewrite our central equations (|6j, Q and d!5l > in 
the near-BEC limit to compare more directly with the case of 
a weakly interacting Bose gas, described by Gross-Pitaevskii 
(GP) theory. It can be shown that (in the absence of FB) 



= A 



4tt v /2^V| H 3 ' 
which, in conjunction with the expansion of Eq. 0, 

m _ /2m\ 3/2 I" J_ A^' 

4irh 2 a s ~ \ h 2 J 8vr [ + 16 /i 2 



(43) 



(44) 
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yields 



a s irnh 2 



2ma; 



(45) 



These equations hold at all T < T c . At T = 0, these equations 
have been shown 1 22, 23] to be equivalent to the results of GP 
theory where one identifies an effective inter-pair scattering 
length a B = 2a s via n B = ilTa ^i / Mb ■ Here n B = n/2 

represents the number density of pairs, [i B = 2/J.+ ti 2 /ma 2 s 
is the "bare" chemical potential of the pairs, and M B w 2m 
the pair mass. 



B. Effects of Feshbach bosons 

We now show that in the presence of Feshbach bosons this 
equation of state is no longer that of GP theory and, more- 
over, there are important differences in the ratio of the bosonic 
to fermionic scattering lengths. As a result of the Bose con- 
densate n® in Eq. (I26> one finds an extra term in the number 
equation 



, (l-Uo/ U*) 2 
' i 
% 



4n^/2rn\|I\h 3 
Combining the gap and number equation yields 



(46) 



in 



4irh 2 c 



2m 
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8tt 



16^ 2 



1()7T 



(l-Uo/V) 2 

53 



(47) 



Solving for u in terms of a* one finds a new equation of state 

H 2 2n 2 g 2 h 2 



/'• 



2ma* 2 m U§ 



tot *4 

-n a c 



(48) 



to lowest order in a*. The second term in the above equation 
derives from the Bose condensate term. The first term i n Eq . 
J46i contributes a term of order a* 7 to this correction 12^1 . 
This behavior should be contrasted with the situation when 
FB are absent, where a B — 2a s . It should also be empha- 
sized that the fermionic scattering length in a model without 
FB is an independent experimental parameter, while here a* 
depends on /i, and must be obtained self-consistently. 



where we have used Eqs. i2l\ and (|20| 
lowest order in na 3 



M* 



2m 1 



After expanding to 



(50) 



From Eq. dl5> we then conclude that, at T c . 

^=^6(O q) T c ). 



(51) 



Equation i51\ reflects the fact that, in the near-BEC limit, 
and at T c , all fermions are constituents of uncondensed pairs. 
From the above equation it follows that (MqT c ) 3 / 2 oc n = 
const, which, in conjunction with Eq. j50t implies 



(52) 



Here Tj? is the transition temperature of the ideal Bose gas 
with Mo = 2m. This downward shift of T c follows the ef- 
fective mass renormalization, much as expected in a Hartree 
treatment of GP theory at T c . Here, however, in contrast to 
GP theory for a homogeneous system with a contact potential 
l2(ill . there is a non- vanishing renormalization of the effective 
mass. 

We turn now to the analysis of the behavior of T c in the 
presence of Feshbach bosons. The (inverse) residue in the T- 
matrix is replaced by 



Z\ — Zn 



[U(2fi g 2 



(53) 



and the h 2 q 2 coefficient B± = l/ (2M*) in fi q is such that 



BqZ 



2M [U(2/i-u)+g 2 ] 



2T2 



1 

4m 



(54) 

Here Zq and Bo are the appropriate counterparts when FB are 
absent. Since Z$ is proportional to the fermionic contribution 
to the density, it is very small in the BEC limit. Using the 
same reasoning as in the previous case, we conclude that the 
ratio T c /B\ is constant with varying coupling. Thus, again, 
T c follows the behavior of the inverse effective mass with, to 
leading order, (T c ° - T c )/T c ° oc a 



*'.> 



IV. CALCULATIONS AT T c 

We turn next to a calculation of T c , which requires that we 
determine fl q [via the T-matrix of Eq. d 1 8 M as a function of 
the scattering length a s . We address the one-channel case first. 
The general expression for 1/Mq in the near BEC limit is 
given by 



21.2 



4:E k fl 2 k 

3m 2 A 2 



(49) 



V. PARTICLE DENSITY PROFILES IN TRAPS 

The differences between the equations of state for the one 
and two channel models will have physical implications in 
the density profiles of particles in a trap. We now introduce 
the harmonic trapping potential V(r) — ^muj 2 r 2 , which is 
treated in the Thomas-Fermi (TF) approximation. In this ap- 
proximation, one replaces u with fj,(r) = fj, — V(r). In con- 
trast to the uniform case, here fi pa ir ( r , T) becomes non-zero 
beyond a critical radius R C (T), where R C (T C ) = 0. In this 
way at T c , only the center of the trap is superfluid, while at 
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FIG. 2: Density profiles at (a) v Q = -1260 and (b) v = -250. 
Solid and dashed lines are at T = with and without Feshbach 
bosons (FB), respectively, and the dotted lines are at T = T c . Inset to 
(a) plots the FB weight in the condensate. Inset to (b) shows the ratio 
of the inter-boson to inter-fermion scattering lengths at trap center. 
Units are chosen such that kp = 1, and Ef = (3N)~shuj = 1. 
Hence the units for Uo, Ef /k%, and for go, Ep/k^/ 2 , are both set 
to unity. For vo — —250 the interaction parameter kFa* s ~ 0.22. 



T = all of the trap contains condensed states. To obtain the 
T = density profile, n(r), we insert /i(r) into Eq. J47i and 
solve for n(r) [Here n(r) refers to the sum of both fermion 
and molecular boson contributions]. The solution is 



n(r) = 16/i 2 (r) 



1 



y/2m\n(r)\a. 
2m 



1 



(1 



Ua_\2 



9l 



(55) 



where we use N — J n(r)d 3 r to self-consistently determine 
fi. Here it should be noted that a* is itself a function of fi(r). 

In Fig. |3 a), we plot n(r) for the parameters Uo — —0.89, 
go = —35, uq = —1260 and N — 10 5 . For this value of vo we 
are somewhat away from the deepest BEC regime as is nec- 
essary to ensure the validity of the TF approximation and the 
fermionic contribution to the density is no longer negligible 
compared to its bosonic counterpart. Indeed, the percentage 
weight of the molecular boson condensate (shown in the inset) 
indicates that by vo ~ —400 the Cooper pair condensate is be- 
ginning to dominate. In this figure we also show the density 
profile as computed in the absence of FB (dashed line), as well 
as the behavior at T = T c (dotted line). Our T = T c curves 
were computed in the absence of FB for both panels. We were 
unable, thus far to find important differences between this dis- 
tribution and that of an ideal Bose gas. Here one has to solve 
self consistently for // pcM >(f)> as well. 



To arrive at a meaningful comparison of the two T = 
cases (with and without Feshbach bosons) we used the same 
value for the effective two-body scattering length a*, obtained 
from the self-consistent calculations of n(r) in the presence of 
FB. The profile without FB bosons is then calculated using the 
familiar TF result near the BEC limit 1 23] at a fixed fermionic 
scattering 1 30] length a*. The same plots are presented with 
vo = —250 in Fig.[2jb), which is further from the BEC limit. 

An important consequence of the one vs two channel prob- 
lems is that for the latter, as the extreme BEC limit is ap- 
proached, the number of fermions is diminished in favor of 
FB. This occurs due to the self consistent adjustment of the 
fermionic chemical potential. Thus the interaction between 
the bosons mediated by the fermions is weaker than what 
might have been expected without consideration of FB. This 
weakened interaction is reflected in Fig.|2]through a compari- 
son between the solid and dashed lines, which shows that the 
trap profiles are narrower in the presence of Feshbach effects. 
Indeed, this could have been anticipated from the above cal- 
culations in the homogeneous case, at T — and T c , which 
deduced a very weak dependence of the bosonic scattering 
length and effective pair mass on the fermionic scattering 
length a*. This finding mirrors the result in a trapped atomic 
Bose gas 1 20] where one sees that a weaker repulsive inter- 
boson interaction leads to a narrowed density profile. Com- 
parison between the lower and upper panels of Fig. [2] shows 
that in both cases (with and withour FB), the profiles become 
narrower as the BEC limit is approached. 

We may use the results in Fig. |2] to obtain a semi- 
quantitative estimate of the bosonic scattering length as, 
based on a phenomenology used in experimental analysis 
{3 S 0- We compare our results to the Thomas-Fermi ap- 
proximated GP equation at T — which yields ng(r) = 

,1B ~Ub^ wnere me inter-boson interaction Ub is connected 

B-. If we fit the 



to the scattering length ag via Ub 
profiles of Fig. |2]to the inverted parabola (C\ — C^r 2 ), we 

may infer that as = ^K?c 2 ' ^ e ra ^° aB / a *s ^ s plotted vs vo 
in the inset to the lower panel of Fig. [2). The same analysis 
applied to the profile without FB yields the familiar 2:1 ratio 
of the bosonic to fermionic scattering lengths. The parame- 
ter as is an important quantity which appears in experiment 
to be considerably less than a factor of 2 times its fermionic 
counterpart. 



VI. CONCLUSIONS 

There are very detailed calculations [31] of four-body 
atomic scattering processes which yield a fixed length ratio 
of 0.6 for the bosonic to fermionic scattering lengths. Pre- 
sumably in the asymptotic BEC regime this constraint should 
be imposed on a more exact theory of the ground state. In- 
deed, a weakness of the mean field approach taken here (based 
on a generalized BCS ground state) is that the "inter-boson" 
scattering length is treated only approximately. This can be 
viewed as related to the inclusion of only one fermion and 
two fermion propagators (i.e, T-matrix) with no higher order 
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terms. 

A strength of this approach, over the counterpart atomic 
calculations is that many body physics associated with the 
broken symmetry or superfiuid state is included. A compar- 
ison of the one and two channel calculations presented here 
shows that, through many body effects, the nature of the con- 
densate enters in an important way to determine the equation 
of state and, thereby, the ratio of the scattering lengths. It 
should also be stressed that the many body context with which 
the number 0.6 is to be compared is different from 2.0, which 
is the number generally assumed in the literature. This dif- 
ference appears once one includes Feshbach bosons. Indeed, 
while it depends upon the magnetic field or detuning, in the 
near-BEC regime this ratio is found to be significantly less 
than 2.0. This is a central point of the present paper. Future 
experiments will be required to determine, whether, as some 
have presumed, the ratio of 0.6 applies to all detunings, or 
whether this number is variable as argued here. 

In this paper we have shown that superfluidity of fermionic 
atoms in the near-BEC limit is in general different from Bose 
superfluidity (as described by GP theory). We have compared 
one and two channel models and find that differences from 
the GP picture for each are associated with the underlying 
fermionic character of the system. This shows up at T ^ for 
the one channel system and at all T for a Hamiltonian which 
includes Feshbach bosons. 

Essential to our approach is the presence (at all T ^ 0) 
of non-condensed pairs which contribute separately to the 
fermionic gap parameter A, in addition to the condensate 
term. Noncondensed pairs or "pseudogap effects" have been 
generally ignored in the literature 1 24] and we have shown 
elsewhere 1 19] that they are negligible only in the strict BCS 
regime. Our inclusion of these effects corresponds 13211 to 
a self-consistent Hartree treatment of pairing fluctuations in 
BCS theory, in which condensation naturally occurs in the 
presence of a finite excitation gap ("pre-formed pairs"). Our 
approach also makes evident the strong analogies between this 
generalization of BCS theory and Bose condensation, where 
both condensed and non-condensed bosons must be prop- 
erly characterized. In the two channel model, non-condensed 
bosons should be viewed as consisting of a strongly hybridized 
mixture of fermion pairs and Feshbach bosons. In the one 
channel model they result exclusively from fermion pairs. 

A key aspect of this work is in the comparison we have 
presented between one and two channel models. Differences 
arise because of the nature of the condensate. A fermionic 
or Cooper pair condensate (A sc ) enters into the gap equation 
while a bosonic condensate (n[|) enters also into the number 
equation. When there is an appreciable fraction of bosonic 
condensate (as in the near-BEC and BEC regimes) these dif- 
ferences will be apparent. A distinction between the one and 
two channel models is relatively unimportant once the bosonic 
condensate is negligible (usually when the scattering length is 
large and positive or alternatively negative). These differences 
are associated with physical properties such as the equations 
of state and the density profiles in a trap. Feshbach bosons 
can, in effect, collapse much more completely to the center of 
a trap than can fermion pairs. This effect can be inferred from 



the ratio of the scattering lengths, and is related to the fact that 
in this limit, fermions are essentially absent, since the number 
equation constraint can be entirely satisfied by populating the 
bosonic state. As a result, these bosons are close to ideal. 

For the relatively broad Feshbach resonances currently un- 
der study in lithium and potassium, differences between the 
one and two channel models are presumably only important 
in the BEC and near-BEC limits studied here, where there is 
an appreciable bosonic condensate. It will be interesting to 
study narrower Feshbach resonances where these differences 
may persist into the unitary scattering regime. 
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APPENDIX A: CONNECTION BETWEEN x(Q) AND THE 
STANDARD GAP EQUATION 

The expansion of the T-matrix at small four-vector Q = 
(f2, q) can be written in the following form: 



t 1 (Q) = Z (fl - £l q + (lpa iT + iTq) 



(Al) 



Therefore, the condition for the divergence of the T-matrix at 
zero Q is equivalent to 



l^pai 







or 



u- 1 + x (o) = o. 



(A2) 



(A3) 



We now show that for a proper choice of the pair susceptibility 



x(Q) 



K 



G(K)G (Q-K) 



(A4) 



Eq. JA2t is equivalent to the BCS-like gap equation of Eq. 
provided also the Greens function is of the BCS form 



G(K) = 



lOJ n 



(A5) 



Just as in BCS theory, we associate this dressed Green's func 
tion with self-energy 

-A 2 G, 



-K) 



(A6) 



We calculate \{Q — (fflm: l)) by performing the appro- 
priate Matsubara sums following standard procedure 10311 . At 
Q = 0, we have 



X(0) = 



5 (^«) 2 



6 



l 



El 



E 1 



2/(£ k ) 



2E V 



(A7) 
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and Eq. ( IA3I becomes 



gives 



^ 1+E l-2/(i? k ) =|| 



2Eu 



(A8) 



E 



1 1 - 2/CgjQ 
2e k 



2£ k 



(A10) 



which, in conjunction with the two-body scattering equation 



Anh 2 a s 



(A9) 



Thus, we have demonstrated that, if we presume Eqs. ( IA4I > 
and ( IA6> , then Eqs. ( IA2t and ( lAlOi are equivalent. 
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